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Introduction 

This  report  Is  e  tunoMry  of  the  results  of  the  reseerch  carried  on 
during  the  year  under  contract  Mo.  Ml40( 70024)707^^.  The  purpose  of  this 
work  was  to  attempt  the  development  of  mathematical  techniques  for  the 
determination  of  pressure  maaima  in  acoustic  transducer  arrays.  In  general 
two  different  methods  of  attack  were  suggested*. 

(1)  the  study  of  the  problem  using  the  techniques  of  generalised 
harmonic  analysis^  moment  generating  functions,  and  expansions  in  terms  of 
special  functions,  and 

(2)  the  development  of  an  analog  technique  for  evaluating  the  integral 
representation  of  the  solution. 

The  report  is  divided  as  follows: 

Part  I  The  techniquea  of  integral  transforms  are  used  to  develop  an 

integral  representation  of  the  velocity  potential  for  pistons 
of  arbitrary  shapes.  Specific  attention  is  given  to  the 
circular  piston  with  the  result  that  the  expression  for  the 
velocity  potential  is  sin^lified.  The  result  is  then  compared 
with  that  obtained  by  Louis  V.  King  as  published  in  the  Canadian 
Journal  of  Kesearch. 

Part  II  The  methods  of  Part  I  are  applied  to  a  rectangular  piston  and 

a  tentative  solution  for  the  velocity  potential  is  obtained 
^ich  is  analogous  to  King's  result  for  the  circular  disk. 

Part  III  A  brief  discussion  of  an  alternative  formulation  of  the  problem 
is  presented  and  an  expression  is  derived  which  relates  the 
velocity  potential  at  a  given  fixed  point  to  the  potential  at 
its  projection  on  the  (xy)  plane. 

Part  IV  The  possibility  of  analog  computations  with  the  maximum  use  of 

standard  analog  equipment  is  investigated.  Here  the  problem 
is  restricted  to  an  arbitrary  array  in  an  infinite  plane  baffle 
and  to  pressures  on  the  surface  of  this  plane.  Included  are 
numerous  graphs  and  pressure  contours  obtained  by  digital  compu¬ 
tation.  These  results  are  compared  with  the  computations  of 
Sherman  and  Kass  as  found  in  U.  S.  L.  Report  Mo.  495  and  will  also 
serve  later  as  a  check  for  the  new  analog  approaches. 
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Part  I 

The  Basic  Equations 

We  atart  with  the  fact  that  the  velocity  potential  resulting  from  any 
source  distribution  on  an  infinite,  rigid  plane  8,  is  given  by  Rayleigjh,  J.W.S., 
The  Thaorv  of  Sound.  Vol.  II,  MacMillan  and  Co.,  London,  19^0,  p.lOT. 


velocity  potential 

S  integration  of  entire  plane  S 

f  distance,  source  point  to  field  point 

boundary  condition,  l.e.,  normal  velocity  distribution 
on  the  plane 

(3)  acoustic  pressure  * 

Green's  Function  for  scalar  Helnholts  equation 
for  the  plane. 

To  avoid  carrying  along  excess  symbolic  baggage  we  shall  discuss  the  nature 
of  y  on  the  plane  S. 

Equation  fi)  becomes 

•*  •  ^ 

vhere 

(5) 

points 

source  points 

-It ,  u»/t  =L 

^  velocity  of  sound  In  the  medium. 

Points  on  the  plane  are  characterised  by  S  •  0. 
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Integral  Trantforms  of  the  Preasure  Distribution 

Confining  ourselves,  for  the  moment,  to  the  velocity  potential  on 
the  plane  ^  -  0,  we  have 

The  first  step  Is  to  take  Its  Fourier  Transform 


By  the  Fourier  Inversion  Theorem  we  heve  an  integral  representation  of  '/'a, 

/ ►  rf  /■  \ 


Applying  this  Idea  to  equation  (4)  gives 


(8) 


i 


cx- 


which  reduces  to 

00 

(9)  ’^'7'  ^  )  wCa.^j 


when  we  define 


g _ 

y[U^lyT{flo^ 


The  evaluation  of  the  integral  defining  ddpends  on  a  theorem  due  to 
Bochner  and  Is  the  subject  of  the  next  section.  For  the  moment  continue 
irlth  more  physical  considerations  and  Introduce 

fi(r,  »l  ;  4  )  -  1  -  V/*  (1, 1;  *  j  \  K- 

for  the  square  of  the  magnitude  of  the  velocity  potential  %ihich  Is  proportional 
to  the  square  of  the  magnitude  of  the  preasure.  Because  our  interest  is  in 
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limitations  in  maximum  driving  power  due  to  the  onset  of  cavitation^  it  is 
only  the  magnitude  of  ^  or  that  is  of  interest.  The  introduction  of 
Fg  is  the  crucial  step  in  the. analysis.  Continuing  with  a  substitution  of 
the  integral  representation  of  we  have  . 

X  *  I  /  *  >  as-;  £CAtX')  J  III/ 

L®  t  ^  ^  *  I  *  *  ^  * 

**  jtYAV*  *  ^ 


We  then  define 


and  are  able  to  write 

Finally  by  substitution 

\  Liy,  *)  -  «,■»  <  *■>«.'  <^*-A.>“  ■-»,  »;  7 ‘  ^w<'«.7 

the  Fourier  transform  of  the  absolute  pressure  in  terns  of  the  velocity 
distribution  and  function. 

A  similar  analysis  will  yield  the  corresponding  expressions  for 

)  «nd 

Evaluation  of  the  Transforms  by  Bochner*s  Theorem 
As  before  define 

( 1 1 )  ^ = /XT  v's  ^  j  f  J  5  • 
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Then  we  have  the  Integral  representation 

Again  application  of  this  idea  to  equation  (4)  gives 

( 13)  Si  (*, A»;  * ) '  - ».■>)  * 


It  then  becomes  convenient  to  define 

e*" 


(lU) 


4=,  dNi-ir 

*75? 


Both  “fc  )  And  •*.  )  may  be  evaluated  by  a 

Theorem  due  to  Bochner,  A>g<>  S.  Bochner:  “Lectures  on  Fourier  Integrals", 
Annals  of  Mathematics  Studies  No.  U2,  Princeton  University  Press,  1959; 
especially^  43,  “Trigonometric  Integrals  in  Several  Variables".  The 
reference  Just  given  is  a  translation  of  S.  Bochner;  “Vorlesungen  iiber 
Fouriersche  Integrale";  Akad.  Verlag.,  Leipaig,  1932,  reprinted  by  Chelsea, 
N.  Y.  1948. 

He  begin  by  stating 

Bochner 's  Theorem  (Ref.  "Fouriersche  Integral",  pp.  186- I87). 

If  the  absolutely  integrable  function  f(x,,...,  ^1,  )  depends  only  on 

the  quantity  T •  *'»*’ **'  i.e., 

U^n  )  -  <P  "  ^  V  ') 

then  the  function 

Ji*,  .  •  •  ■  i  . /  f  ^ 

in  a  similar  manner  depends  only  on  the  quantity. 

Furthermore,  one  can  replace  the  original  -fold  integral  for  j(®t  )  by  the 
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Simpler  expression 


where  (t)  Is  the  Bessel  function  of  the  /»>  order.  If  one  further  introduces 

L  ^ 

the  quantity  s  *“•'  *  one  discovers  that 

1. )  For 

2. )  For  *.«,/  J(*).(-.rTr“4""  •'Vj.-i 


From  eq.  (lU)  we  have 


■*)  *  '  '/arrlji* 


.alii. 


JtJtiS 


To  apply  Bochner's  Theorem  to  this  we  make  the  following  substitutions: 
JVA,  ot,  f'XV-k, 

VF>  «  Vj.  t  ^  *a 

jwp  •  3  •  *» 

¥*i  ^ 

By  Bochner's  Theorem, 


where 


r  a  k,*-  #•  f  X,* 


Igs  i  9  :.  yn  ai  t 


JU)a  -  fir  -»/d5,  J  casC»^y>)  JjD. 

ixpression  becomes 


-5- 


PaBU  MATBBMATIOAL  LaBOBATOBOM,  lNOOBr<»ATBD 
BBDFOBO  BOAD  *  CABUILB,  MAIIACnTOMm 


where 


1  A  * 
•!. 


The  problem  that  remains  la  to  evaluate 

Let 

J/j.,  I 

The  Integral  expression  may  be  %nrltten  In  the  following  form: 

O  * 


From  Oberhettlnger,  Fritz:  “Tabellen  cur  Fourier  Transformation", 
Sprlnger-Verlag,  Berlin  1957>  p.  20  • 

W4*/. 

'-Ikl-  V^-4-1'-  ’Vvl] ^ h 

c« «4/.)  J>.  =  -  I  (»?/.)  Jj>  = 


/  r** 

y*is  I 


From  the  Vorwort  of  the  preceding  reference,  the  function 

U)  s;»,  (KyJ 

after  a  Fourler-Slne  Transformation  is 

f<x)«**/7rj^  »/> (BJf  J 
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Then  r  r  .  i  • 

r  f  f(x'^dx'  .  A  2  6. 

*^0  ^ 

By  definition  ^ 

Applying  this  to  our  Integral 


Therefore 


“  '/af^  ^  ait9(%jg)  si*  (f^j>)  )■ 

>re 

T(-fc  s  )  - — -■  ■>  -f 


Equation  (l6)  becomes 


-srl  ('Ah 


e  2 


_ f  iLi- 

L'fe'-s.  J 


uhere  ^  ^  3- 

Returning  now  to  equation  (lO)  and  applying  Bochner's  theorem  using 
the  same  substitutions  as  before 

(18)  ,.y,*)  .  0-.  (^r./.)  Jj.. 

From  Oberhettlnger,  Fritz;  "Tabellen  zur  Fourier  Transformation”,  Springer- 
Verlag,  Berlin,  1957,  P-  6l  _  1.1-1 


\  eoiCijs)  Jf  * 

Also  from  p.  16^ 

5  = 


Therefore 


r^a-*‘ 


ox"*  ^  /3" 

*>/3I 

«  Z. 

ie>V^ 


4%  J* 
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where 


X,\ 

s  =  c  (  X  -#->*  y  • 


Compatabillty  Check  between  K„  and 

and  are  related  and  this  relation  furnishes  a  valuable  check 
on  some  rather  tricky  integration  techniques  used  to  obtain  them.  The 
relation  between  them  is 

-«o 

To  prove  this,  begin  with 

^  s. (A - 

Since  the  exp  [-AvL(\x*/i^f)  ]  factor  does  not  contain  let 


L.  *  - 


fU*' 


TTfc 


S,  (-fe  ^  =  (^tr;  V 


Therefore  the  relation  that  must  be  shown  is 

U  =  •^■‘4  J*- 

The  right  side  of  the  equation  is  equal  to 

_£fV-ca7r)‘av/.Vjrv  “  ^  ^ 

An  analysis  of  the  first  integral  yields 

r**  ^  _  J  f*  *** 

*  J  •*=  ‘  -  'a  • 
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Let  SubiCltution  gives 

'fr  ■ 


“•B 

Let  */a  ®  ^  Therefore 

j/  .  ■  ‘L  r ju  • 

From  Ryshlk^  1.  M.  and  Gradsteln,  I.  S.>  "Sumnen;  Produkt  und  Integral 
Tafeln**,  Veb.  Deutscher  Verlag  der  Wlssenschaften,  Berlin,  1957*  No.  1. 513*2 

/ 

Therefore  the  first  Integral  has  the  value  sero. 


L,w  A-4±I- 


* 


=  o. 


Since  the  second  Integral  Involves  an  even  function 


With  the  following  substitutions: 


.  xJx  «Cai»)N  Jv, 


^aT>V 


[.if  S.(4-«)t-- 

*  *  y  x‘-Ca-ir)‘aV/iV 


In  the  previous  section  the  function  was  defined  as 

j  (Hi-  x)=  Vt-  J 

*  * 

We  shall  show  later  In  the  Appendix  that  can  be  evaluated  by  the 
direct  substitution  of  this  quantity  Into  the  Integral.  Here  we  shall 
treat  Sg  as  the  Dirac  delta  function  and  because  of  the  limits  Involved 

o  “fe  « 


,  ,=s==  -6>irr 


y, 
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The  Integral  3  then,  has  the  value  •; - -  ■  -  or  -  * 

when  w  and  thus  la  equal  to  Lg* 

Comparison  with  Klng*a  Reaulta  for  a  Circular  Platon 

Louis  V.  King  wrote  an  excellent  paper  "On  the  Acoustic  Radiation 
Field  of  the  Pleso>E lac trie  Oscillator  and  the  Bffect  of  Viscosity  on 
Transmission",  Canadian  Jour,  of  lasaarch  II,  (1^93^)*  He  obtains 

the  expression: 

ai  J 


for  the  velocity  potential  due  to  the  disc  of  radius  "a"  oscillating  with 
slaqile  hamonlc  notion  and  naxi—  velocity  x.  Ihasa  synbols  need 
translating  for  comparison  with  our  work  and  this  will  be  done  after 
obtaining  the  sane  result  by  our  nethod. 

For  the  King  case 


W 


W. 

o 


otWi>wise 


Then 


But 


- J 4)  Jx  J.J  . 
\  =  c  ^ is 


%)here  L^ 


A 

V-J, 


Since  does  not  contain  x  and  y  then  by  substitution 


*  Bead  [  P;  condition  on  P]  "Che  set  of  points  F  an^  that  the 

condition  on  f  is  sstisfIniP. 
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This  expression  has  the  correct  symmetry  for  an  application  of  Bochner's 
Theorem  where 

With  the  proper  substitutions  and  application  of  the  theorem 

CA,/.,!);  .  W.  [«T  I  3-.  Jj,] 

The  expression 

^ j  ^ /)  J/»  -  •* 5  ^  ‘=!p 

which  Is  evaluated  by  No.  335*1  from  "Tables  of  Integrals  and  Other  Mathe 
matlcn^  Data"^  Dwlglit^  Herbert  Bristol,  MacMillan  Co.  New  York  to  be 


f 


I  ^  ^  J  ^  A, 

a  stli«yuise- 


Thus 


A ' 


W.U.  ^  3-,C«;a). 

From  Equation  (12) 


Substitution  yields 
oe> 


i»r<- 


C»?v'T 


J-tf 


where 


(aurCAS/*") 


The  first  step  in  the  evaluation  of  q-, ^ H,  5;  4 )  is  to  consider  the 


Integral 


J 


A 

r 


TTc 


-i* 


\ri; 


6,  (“fe  )c 
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Since  S,»  ^  a.  (a-^)  y  have 

f‘->  - TT  f  -r— -"'jf 

J  -ie  -S^~  (»v)  j)*-  y  y/^  t  * 


jc 


where  we  will  denote  the  first  Integral  by  A  and  the  second  by  B. 


-CaTr)\ 


Jif 


_  ^  fl_ — :? -  cos  <iT>. 


Let  A  « 

x‘=  (■ar;‘v‘- 

<1>1  =  Jtv  Jt< 

d*.  ib*^. 

Then 

From  "Tabellen  zur  Fourier  Transformation" ,  P*3 

•i^ 

x‘;''coaCa^J  dx  -  WsA  sm  (a;;. 
Then  ^  _  sin  when  $  P  O. 

V**-  s. 


Let 


_f - Ei: -  (■€ #•  (•»)•->** 

y^TuiFirir 


- 


_  ^  f-  5J-4  -\/5,  Af4r;v- — ^ 


Am 

dx  s  (*ir)*'o>  Jv 


xix  -  Ji) 


V  =ijIE5 
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Then 


L  rrrrr. 


Treating  Sac'll- x)  «■  we  did  In  the  previous  section 

■ft  >  \rir 

I  \/  ■*“_ 

B  a  L 


c*x  SIV-^  j, 

/ipr^ — 

0  -ft  -c. 

The  Integral  under  consideration  then  Is 


A  +  B  "  ^*** _ ^  ~  ^  Cas  5 \f‘fc  -  g».  -ft  > 

\rvTT^ —  >nFnz —  ' 


'/S‘-  5* 

Therefore 

•0 

a>B«  ^ 

a;«r.*) 

Again  *(»»)*■ 

and  we 

^>rsZ  . 


rrt. 


t'At 


which  Is  suitable  for  the  application  of  Bochner's  Theorem. 


Let 


X  a  ai»l 


dX 


d», 

*-r 


I.  . 

=  — » 


and  we  obtain  by  substitution 


Application  of  Bochner's  Theorem  yields 
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If  we  let  ^  ^  -  / r‘-  ■*** 

^  H ,  -fc)  -  a  VJ.  j  j, («*»)  T,  (fj.)  Jr 

wliich  is  comparable  to  King's  results. 


Part  II 


The  next  case  to  be  considered  Is  the  investigation  of  the  velocity 
potential  for  a  rectangular  piston,  centered  at  the  origin  with  half-sides 
"a",  "b"  as  shown  in  the  figure. 


We  start  again  with  the  equations  derived  in  the  previous  section: 
(1) 

(2) 


For  he  pressure  boundary  condition  v:'  take 

r  1  for  points  (x,  y)  inside  the  rectangle 
•  I  0  for  points  (x,  y)  outside  the  rectangle. 

Substitution  in  equation  (l)  yields 


rr^ 
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After  integration  the  expression  becomes 


.  /.  .  -  \  r  ^  ^  5i'>i  (••sm'*)  S'lViCsirAa) 


The  general  integral  representation  for  the  velocity  potential  was 
found  to  be 

f  *•>•1  ) 

•  •• 

Substitution  of  equation  (5)  into  this  expression  gives  us  the  velocity 
potential  for  the  rectangle 


(7)  s'"  JU,. . 


-•* 


«Vi  VS 


The  evaluation  of  the  second  integral  was  carried  out  in  detail  in  Part  I 
and  found  to  be  equal  to  the  following  expression  when  >  anrr  V  A*-  •*■/** 


(•■*)*■(  A' ^-M*) 
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After  substitution  and  separation  of  variables  our  expression  Is  of 
the  form 


S;i,C«Tra.A) 

-n“A 


-16- 


PaMDI  MaTUMATWAL  LaMBATOBIII,  iNOOVOAAnD 
MWOID  aOAD  •  OABLULB,  MAtfAOBUtBTTC 


or 


J/i, 


From  Campbell  and  Fpater^  Fourier  Intagrala  for  Practical  AppllcAttona . 
19^  printing,  p.  Ill,  pair  Mo.  868 


S(,). '/r 


uhere  end  d'  ere  complex  numbera,  not  Infinite  and  the  real 

parta  are  not  leaa  than  0. 

Applying  this  to  our  expreealon  ^ at.  ~  f  ^  Itt^ 

■‘L„  =  ■»  K.  [i  !/**•_  ]  . 

Therefore 

A.  aJ  K.  J  J":- 

-•o 


similarly  let 


Introduce 

Then 


B  *  g  •  ^ 

_i,  V  a*  —  Cjotm)** 

a*”* 'fc*'—  . 

-i  A  iT**-  ^ 


Applying  the  same  method  aa  before 

B  »  A  J  /€*♦  5*-  J  • 

The  expression  X  *  '/a  A  'a  'B  then,  becomes 

— K+'b  ^  *1“  b 

I  =  J  -  laiTA)*"  5  4c  —  J  K, t*  t/c* 

Therefore 

X«J  K^lW||«--UirA)‘>)c»A|*]  dc  . 

\-V 


BSDM»D  mo/a  •  tUkmUBIM,  MAMAOBUtSm 


Substitution  into  aquation  (8)  yialds 


SiM  (air«,>^ 


iJlf 


J  K.DW-CairX)^  4c. 


Our  final  axprassion,  than,  for  tha  valocity  potantial  of  tha  ractangle 
is  givan  by 

••  afi 

(io)v,(SA,S;i)- J  (  ‘  c*»x)»-Ve»^5*  JclciA, 


This  is  only  a  tantativa  solution  to  tha  problaa  of  tha  rectangle  which 
is  analogous  to  King's  solution  for  tha  clrela.  Further  study  is  necessary 
in  order  to  dataruina  whether  this  double  integral  can  be  reduced  to  a 
sinplar  foru. 


Part  III 

In  this  section  wa  discuss  briefly  a  certain  alternative  formulation 
of  tha  problaa.  Work  on  this  has  not  yet  progressed  to  the  point  where  it 
is  possible  to  make  a  definite  stataawnt  as  to  its  value,  but  we  feel  that 
it  Bwrits  further  consideration. 

Let  D  be  a  doewin  In  the  plane.  The  basic  expression  for  the 

velocity  potential  for  the  domain  H  Is,  from  Part  I, 

(1)  *“  />• 


In  which  la  the  boundary  pressure  distribution  normal  to 

t. n)*  ♦  . 


By  a  change  of  variables,  we  may  take  ^  m  & 
the  fora: 


,7-0 

4x4^ 


so  that  (1)  assumes 
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.  Taking,  as  before,  wCn.’j)  •  < 


in  which  X  «  >J  X*''#* 
the  point  is  in  the  domain  P  and  W(x,  ^)»0  otherwise,  we 


if 


have 


(3) 


where  now  the  double  integration  is  extended  over  the  domain 
_ _ 

Let  />a  H  X*"  ¥■  y*"  .  Then  the  function  ^  is  constant  on 

circles  of  radius  Ji  centered  at  the  origin.  This  means  that  the  integral 
(3)  My  be  evaluated  by  summing  up  the  contributions  to  it  from  annular  bands, 
centered  at  the  origin,  in  which 


/**  is  approxisMtely  constant.  Thus, 


referring  to  the  figure,  we  have  that  the  contribution  due  to  the  annular 

band  ABCD  is: 


at  the  distance 

from  the  origin,  lie  define  the  function  to  be  that  fraction  of  the 

circumference  of  the  circle  of  radius  |0,  centered  at  the  origin,  which  is 
intercepted  by  the  domain.  It  Is  easy  to  see  that  ,  to  that 

the  expression  for  the  contribution  due  to  ABCD  becosms 

^ABCP  r  » 


The  values  are  zero  for  those  p  such  that  the  circle  of  radius  p 

does  not  neet  the  domain  1>  .  Summing  up  the  contributions,  we  have,  using 
the  fact  that  ^  Jp  m  r  Jp 


-19- 


Pams  MAtnuiATioAL  LAMNUTOUst,  Imoobmbatbd 

BSWOBS  BOAB  •  CASUWJI.  MiASACISOTlTTt 


and  then^  since  ^  f  -  <  *“ 

(4)  *  -  J 

f 


cir 


Tills  is  the  basic  expression  In  our  new  foruiulatlon.  We  see  that  the 
problem  of  evaluating  the  velocity  potential  for  a  given  domain  splits 

Into  two  parts: 

I)  find  the  function  for  the  given  plane  domain  Tills  Is 

the  purely  geometric  part  of  the  problem; 

II )  carry  out  the  Integration  In  (4).  This  is  the  analytic  side  of 
the  problem,  and  the  more  difficult.  Without  going  Into  too  much  detail,  we 
give  here  a  few  brief  remarks  concerning  each  of  these  two  parts. 

I).  If  ^  Is  a  disc  of  radius  a,  centered  at  the  origin.  It  Is  easy  to 
see  that  ■5'^(  p3  *  for  6*^—  ,  end  •  o  otherwise.  In  this 

case,  the  Integral  (4)  becomes: 

Jr 


and  this  may  be  explicitly  evaluated  to  yield: 

4^C'5;‘fc)=  'film  )  ^ 


-ciS 


-  c. 


This  result  agrees  with  that  of  S tense l-Brosse,  Leltfaden  sur  Berechnun:.. 
von  Schallvorgingen.  Springer-Verleg,  19^8;  Part  II,  4a,  p.  75 • 

If  D  is  a  rectangle,  with  '*half-sldes''  k  >  O' Z.  b,  as  in  the  figure, 

It  may  be  shown  that  ’$‘^(^3  1>  given 

by  the  following  expression: 


1 
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Again,  if  P  is  an  annular  band,  situated  as  shotm,  of  width  J  , 
central  angle  A  and  inner  radius  A  then 


KSf)  -  A 


o 

>u 

Ai-  f>  A  *"J 

0 

f>  > 

For  certain  other  domains  P  of  simple  shape,  the  function 
can  be  given  explicitly.  Of  course,  it  goes  without  saying  that  if  D  is 
irregular  it  will  in  general  be  impossible  to  give  an  explicit  formula 
for  .  In  this  connection,  it  is  worthwhile  to  mention  that  the 

function  !■  additive,  i.e.  if  1>,  ,  are  two  disjoint  domains,  and 

their  union,  then  =  ^i> C/»)  ^ J*)  '  I" 

then,  a  given  domain  P  of  irregular  shape  miy  be  approximated  by  the  union 
of  non-overlapping  domains,  P,,^>  each  of  which  is  of  a  simple 

regular  shape.  Then  (ye)  ^  ***  ^  is  an  approximation  to 

Remark:  Another  interesting  property  of  is  the  following:  If  A(p) 

denotes  the  area  of  the  domain  'P  then 


»  Air  J 
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II ).  If  D  it  A  circle  centered  at  the  origin,  we  have  seen  that  the 
Integral  may  be  explicitly  evaluated.  In  the  case  of  the  annular 

domain,  alto.  It  It  possible  to  carry  out  the  Integration  In  (4),  using  the 
expression  given  above  for  However,  If  b  Is  a  rectangle,  we  already 

run  up  against  a  rather  formidable  analytic  problem. 

If  we  substitute  In  (4)  the  expression  for  given  above  for  the 

rectangle,  and  simplify,  we  obtain  expressions  Involving  Integrals  of  the 
form 


(5) 


I 


/  1*^ 


e  Jr 


yf  Jf'- 

In  which  c  >  A  ^  A. 

This  type  of  Integral  seems  very  difficult  to  evaluate.  One  potentially 
successful  method  of  handling  It  may  be  the  following. 

Let  denote  the  Integral 

r*  I 

where  f  Is  an  arbitrary  function,  subject  only  to  the  condition  that 
^  f(x)  Jx  exists.  If  ve  set  V  t*-  J  •>  we  obtain: 


(6) 


ffl) 


From  Campbell  and  Foster,  as  above.  Part  II,  we  have: 


(7) 


Substituting  this  Into  the  expression  (6),  and  formally  inverting  the  two 
Integrals  gives: 


(8) 


/i^ 


sK«)c’*^  J*. 
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Now  ^  Js  «  -  V*  9)  as  Is  evident  from  the 

definition  of  1.  Thus  (8)  becomes: 

■o 

(9)  Jk.  1(0, ,)  J  J  • 


Integrating  by  parts  we  obtain:  ^ 

•A 

From  Higher  Transcendental  Functions.  Bateman  Hanuscript  Project, 
Erdelyi  et  al;  Vol.  II,  p.  79,  we  have: 

‘'/j,  K.t*)  -  -  K /»)■ 

Using  this  in  (lO)  yields: 

•o 

(“)  1  ^  [k.  (;  - 
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Since  (Hi^er  Transc.  FuncClona,  ibid.,  p.  23)  V^«^**)*  0,and  alnce 
X(o-^^  remains  bounded  as  we  see  that 

\.K<.-i'rY^)Xo>,f)]  =  o 


and  (11)  becomes: 


This  last  equation  is  our  desired  result.  Its  slgnlglcance  Is  that 
It  relates  the  evaluation  of  X  it)  "Ik)  to  that  of 

If  for  fC»)  we  take  for  a  given  domain  then  X^S,  ♦) 

Is  just  and  (12)  specialises  to: 

(13) 


S  f 


Thus  the  velocity  potential  at  the  field  point  '^)  with 

propagation  constant  '&  may  be  expressed,  by  this  last  relation,  in  terms  of 
the  potential  for  t  ^o. 

If  for  K'r)  we  take  sin"'  ^/v*  ,  as  in  (^),  we  aae  that  the 

evaluation  of  (3)  reduces  to  the  evaluation  of  the  Integral 

B 

f  -I  - 

J  ^ 

A 
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which  is  the  form  (5)  takes  for  •  o.  This  latter  may  be  transformed 
(we  do  not  give  the  details  here)  into  an  integral  of  the  form 

5  i'C^)  where  f  is  the  Fourier  transform  of  sin  and  H 

is  the  kernel  function  for  an  Integral  transform  of  the  type  of  the  Bessel 
transform. 

Evidently^  these  formal  manipulations  need  to  be  supplemented  by  much 
more  careful  analysis  in  order  to  reach  a  satisfactory  end-result. 


Part  IV 

The  Basic  Equations 


The 
velocity 
plane  S: 
(1) 


starting  point  is  again  the  basic  integral  representation  of  the 
potential  y  due  to  a  source  distribution  on  an  infinite  rigid 

-;*r 


vp  s  ... 


s 


AS 


y  .*.  velocity  potential 

5  .>.  integration  extends  over  the  entire  plane 
V  .».  distance,  source  point  to  field  point 


(2)  W  s  ^  normal  velocity  distribution  on  the  plane. 

(3)  acoustic  pressure. 


— - —  Green's  Function  for  scalar  Helmholtz  equation, 

'  for  the  plane. 


For  a  sound  radiation  field  having  a  single  frequency  component  'P 
is  proportional  to  acoustic  pressure.  It  will  be  convenient  to  deal  with 
the  velocity  potential.  This  is  a  function  of  position  P  of 
the  field  point  and  the  propagation  constant  “Ic  =  ••"A  .  He  thus  avoid 
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carrying  along  aome  constants  that  add  nothing  to  the  discussion  at  hand. 
If  1>  the  velocity  at  every  point  In  the  plane  of  a  baffle  and 

pistons  and  If  (  )  are  the  coordinates  of  an  arbitrary  field  point  'P 

on  the  plane  then 

(4)  )j  WU.^)  — p-  Jx 

where 


(5)  ^  ^  • 

Here  we  take  S  ,  considering  only  points  on  the  plane  T-s  o.  In  this 
special  case  It  Is  convenient  to  Introduce  polar  coordinates 

^  •  r  Snb 

and  consider  y  at  the  origin  ^  ■  c>  ■  We  may  then  %nrlte, 

(6)  y  (o,6j  4)a: J  jw(r,b)e*’  Jy  - 

•  0 

Here,  we  have  expressed  the  piston  velocity  distribution  d)  In 

polar  coordinates  with  the  field  point  of  Interest  as  origin.  This 
entails  no  real  loss  In  generality. 

If  we  wish  to  get  the  position  of  the  field  point  back  Into  the 
picture  we  awy  w^' 

air  •• 

(7) 


where, 


V(V,9,  velocity  distribution  in  polar  coordinates 

with  the  field  point  (  i(  )  as  origin. 


However,  for  the  anaent  we  do  not  carry  the  excess  notational  baggage  of 
equation  (7)  and  use  aquation  (6)  instead,  rinallj  let  ee  ayllt  e^aetion 
(6)  Into  two  equations: 
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iT 

(8)  V(f)-  Jw(»-,«)as 

and 

(9)  j  V(r)e.  Jr  , 

# 

These  equations  are  basic  to  what  follows.  In  general  W(y,^)  Is  complex 
and  the  modules  of  w  Is  the  magnitude  of  the  piston  velocity.  The  argument 
of  w  is  the  phase. 

Some  General  Considerations. 

Before  continuing  it  is  important  to  observe  that  the  area  of  any 
single  baffled  piston  may  be  idealized  to  be  an  area  over  which  w(r,&) 
is  a  constant.  Thus  one  may  write 

N 

(10)  w(r.©)  -  Z  4;  «<; 

ia.! 

where 


<Li  .».  complex  number  giving  phase  and  velocity 
amplitude  of  piston  * 

characteristic  function  of  piston  * 

1  on  piston  A 
0  everywhere  else. 


C.Cr, ») 


This  leads  to  the  consideration  of  the  triple  of  eq^atlon8 


(11)  ''*(0  =  ^  j  J® 


(12)  V'iC‘fc)=  -  I  V.(r)e  J 
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and 

N 

(13) 

iai 

The  study  of  an  analog  technique  can  now  be  broken  down  Into  the 
problem  of  generating  and  calculating  each  of  the  expressions 

finally  the  modulus  of  which  Is  propor> 

tlonal  to  the  magnitude  of  the  sound  pressure. 

A  Function  Generator  for 

The  simple  form  of  the  function  permits  It  to  be  generated 

by  an  optical  device  as  follows.  In  front  of  a  cathode  ray  tube  screen  Is 
placed  a  mask  which  Is  scaled  to  represent  the  plane  baffle  In  the  neighbor¬ 
hood  of  the  array.  The  mask  Is  transparent  In  the  region  corresponding  to 
piston  L  and  Is  opaque  everywhere  else. 


Pabkb  Mathematical  Laboratoiubs,  Incorporated 
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The  null-deflection  point  of  the  CRT  (cathode  ray  tube)  is  set  at  a 
point  on  the  screen  which  corresponds  to  the  field  point  )  on  the  mask. 

The  blip  of  the  CRT  is  set  to  scan  in  concentric  circles  about  this  point. 

If  the  blip  on  the  CRT  screen  is  started  at  the  null  point  and  circles 
about  it  in  circles  of  gradually  increasing  radius  the  entire  mask  will  be 
covered  in  a  relatively  "dense"  fashion.  See  Figure  1. 


Fi^ave  A. 


Light  passing  through  the  transparent  region  of  the  mask  is  focused 
by  a  lens  on  a  photo  multiplier  tube.  The  signal  from  the  photo  multipler 
tube  will  be  a  parametric  representation  of  as  follows: 

(14)  d  =  i-f-Pt, 

where  V  .>.  frequency  of  circular  scan 

(15)  r  *  A'f  >>t 

where  is  the  increment  of  radius  with  earl 
scan  cycle.  ZiV  4.^/. 

Accuracy  of  the  Generator. 

In  order  to  measure  the  accuracy  of  the  proposed  c^(r,9)  generator 
a  manual  calculation  was  made  following  the  method  of  the  analog  generator 
for  a  nine  piston  array  for  which  the  sound  pressure  magnitude  is  available 
in  the  Sherman  and  Kass  report.  (Sherman,  Charles  H.  and  Kass,  Donald  F.^ 
"Near  Field  Sound  Pressure  of  Arrays  of  Pistons".  U.  S.  Navy  Underwater 
Sound  Laboratory,  New  London,  Connecticut,  I96I.) 


-29- 


Paiuu  llAnBMATKUL  Labobavomm,  Inoobmbatbb 

BBSVOBD  BOAB  •  OABUtIJI»  ICAtlAOBUtBTTB 


CalculAtionB  of  were  audo  for  thlrtBon  flold  point*  or  caIcu- 

latlon  stations  located  along  the  axis  of  cantsrs  of  the  mlddla  pistons. 

(See  Figure  20).  Reference  circles  were  constructed  corresponding  to  the 
circular  scan  of  the  CRT.  The  IncresMnt  of  radius  vas  chosen  so  that  the 
array  was  covered  In  73  circles. 

Since  all  the  pistons  were  taken  to  be  vibrating  at  the  saae  frequency, 
phase,  and  aBq>lltude  It  was  unnecessary  to  calculate  for  each 

piston  separately.  All  pistons  were  considered  at  once  giving  Ihe 

symmetry  of  the  array  permitted  the  consideration  of  only  the  pistons  (and 
halves  of  pis  tons  )aicne  side  of  the  axis  of  field  points. 

From  the  values  of  thus  determined  the  sound  pressure  amplitude 

was  calculated  by  nuaierlcal  methods. 

Figures  (7)  to  (19)  show  V(r)  in  graphical  form  for  each  calculation 
station  Vi  .  Figure  (20)  gives  the  comparison  of  the  two  Mthods  of 
calculation  In  the  form  of  a  graph  of  the  sound  pressure  aaqplltude  vs.  length 
along  the  axis  of  field  points.  The  array  Is  characterised  by  'fc*' 
using  the  notation  of  Sherman  and  Kaas,  where  au  Is  the  piston  radius. 

It  may  be  seen  that  the  discrepancy  is  of  the  order  of  1^.  Clearly 
higher  accuracy  will  be  obtained  by  using  a  smaller  increment  of  radius  and 
more  field  points. 

Cgggutatlon  Device  for  i|p,C4)  andl> 


Fi|«rc  3 
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Figure  3  shows  the  suggested  device  for  computing 
that,  ViC't)  and  I* 

A.  Scanning  signal 
^oto  I  pv  *lf 


_  Actual  output  with  afterglow  effects  for  high  speed  operation. 

_  Idealized  output  with  no  afterglow  effects. 

B.  Photo  multiplier  signal 

C.  The  switching  filter  gives  a  rectangular  pulse  output  which  closely 
approximates  the  Ideal  photo  multiplier  output  (with  no  afterglow 
effects)  for  a  high  cut-off  threshold. 

D.  The  Integrator  D,  which  is  set  to  dump  once  each  period,  gives  . 

Taking  r  to  be  a  constant  over  one  period,  this  Is  mathematically 
equivalent  to  Xnr»  j  c^(»^,a)dd  by  equation  (lU).  This  is  within 

a  multiplicative  constant  of  v^Cr)  .  Thus  the  output  of  the  integrator 
is  a  sequence  of  pulses  each  one  of  which  starts  from  zero  and  climbs 
to  a  maximum  value  proportional  to  .  See  Figure  5. 


S*«iVc.ViY»^  F'iVtw 


Ofcl-pWt  I 
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E.  The  low  pass  filter  gives  the  envelope  of  the  pulse  output  of  the 
Integrator,  l.e.  a  continuous  reading  of  Vklv)  .  gee  Figure  6. 


•f 


Figure  (» 

In  Figure  6  the  time  scale  Is  over  many  periods  so  that  t  may  not 
be  considered  constant  but  Is  related  to  by  equation  (l^)* 

F  and  F^.  The  oscillators  F  and  F^  give  sinusoidal  outputs  as  a  function 
of  v*  .  The  phase  of  F  Is  ahead  of  the  phase  of  F^  by  "^/a. 


G  and  G^.  The  multipliers  G  and  G^  give  VfC*-)CM  'fcv  and  -YtCOSut  4»' 
respectively.  ^  ^ 

H  and  The  Integrators  H  and  give  i  V.(r)C©s  •&»>  and  - 

respectively.  This  Is  equivalent  to  *7^^^  ^  V^lr)  C»s  *4**  Bnd 

~  j  5i»  "tr  Jr  It  will  be  noted  that  for  V  >  Yt^mn  *  ®- 


Thus  the  upper  limit  of  Integration  Is  equivalent  to  *0  after  r  reaches 
Yw,»k  ’  outputs  of  Che  Integrators  are  signals  whose  final 

values  are  proportional  to  ^  ’^»(V}and  X»i  respectively .  These 

quantities  are  stored  In  a  recording  device. 
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Automatizing  the  Computer. 

The  system  outlined  above  gives  the  real  and  Imaginary  parts  of 
evaluated  at  one  field  point  corresponding  to  the  null-deflec tlon  point  of 
the  CRT  on  the  mask.  The  CRT  scanning  control  may  be  programned  to  cease 
the  circular  scan  after  the  radius  has  Increased  to  a  value  greater  than  > 

a 

relocate  the  null-deflection  point  to  some  point  neighboring  the  first,  and 
start  the  circular  scan  at  Y*  >  0  again.  The  null  point  advancing  device  can 
be  set  to  cover  the  entire  mask  with  a  lattice  of  field  points. 

The  final  stage  Integrators  (H  and  In  Figure  3)  must  be  set  to  dump 
at  the  completion  of  computation  for  each  field  point.  The  output  of  these 
Integrators  will  then  be  a  sequence  of  pulses  whose  maximum  (  >  cut-off) 
value  corresponds  to  and  Jm  for ‘each  field  point 

(r.n)  • 

It  la  clear  that  the  recording  of  takes  place*  several  orders  of 
magnitude  more  slowly  than  the  circular  sweep  of  the  CRT.  Thus  the  question 
arises:  How  long  will  thie  complete  coverage  of  the  mask  take? 

This  problem  has  not  been  studied  In  detail  but  a  brief  Inspection 
suggests  the  following  estimate.  A  standard  short  persistance  oscilloscope 
coating  retains  10^  of  Its  Initial  glow  2  milliseconds  after  a  sharp  cutoff 
of  the  beam,  and  1^  after  -3  milliseconds.  (Czeck,  J.  The  Cathode  Ray 
Oscilloscope.  Interscience  Publishers,  Inc.,  H.  Y.,  19^7,  p.  21.)  If  the 
filter  which  minimizes  afterglow  effects  Is  %rell  designed  It  should  be  possible 
to  operate  the  circular  scan  at  up  to  200  cps.  About  200  cycles  should 
suffice  for  an  accurate  coverage  of  the  piston  region  so  that  one  field  point 
can  be  calculated  per  second.  A  lattice  of  23  x  23  field  points  could  then 
be  covered  In  less  than  eleven  minutes.  Recent  Improvements  In  shortening 
afterglow  time  may  permit  this  figure  to  be  reduced. 

Systems  for  Computing  the  Sound  Pressure  Amplitude  from  'T,:  andli  Data. 

There  are  a  large  number  of  ways  of  computing  the  sound  pressure  ampli¬ 
tude,  I  pi ;  from  the  values  of  ^  and  Xst  •K  •  Many  of  these  are 
now  under  study  for  feasibility  and  practicality.  A  few  of  the 
and  Xw  recording  devices  and  computation  devices  are  outlined  below. 
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1.  Recording  and  Xh»  V'.;  as  distinct  quantities  for  each  field  point 

with  digital  calculation  of  Ipi- 

2.  Recording  ‘^c  and  Xm  dv  as  continuous  level  curves  corresponding  to 
lines  of  field  points  with  graphical  or  digital  calculation  of 

3.  Recording  i?*  fi  and  Vl;  as  continuous  darkening  of  a  photographic 
film  corresponding  to  the  CRT  mask,  l.e.  corresponding  to  the  region  of 
the  array.  The  light  source  which  exposes  the  film  could  be  restricted 
so  as  to  expose  a  small  region  of  the  film  corresponding  to  the  field 
point,  and  to  advance  as  the  field  point  does.  The  darl.enlng  of  the  fllu 
would  then  represent  the  *4:  or  the  X%  field  In  the  neighborhood 

of  the  array.  Optical  devices  using  these  films  could  perform  the 
multiplication  and  suiBBBtion  In  equation  (13)  for  the  real  and  imaginary 
parts  separately  and  then  give  the  modulus  1*^/  by  vector  addition  of 
cross-polarised  components.  Thus  Ift  can  be  obtained  by  entirely  analog 
means.  The  optical  system,  however,  diminishes  rapidly  In  accuracy  as 
the  number  of  sets  of  pistons  with  the  same  amplitude,  frequency  and  phase 
increase.  The  optical  system  should  prove  to  be  useful  for  accurate  data 
on  simple  arrays  and  qualitative  data  on  complex  arrays.  The  latter 
should  be  of  sufficient  accuracy  to  identify  and  locate  pressure  peaks. 

Conclusions  and  Recommendations 

This  report  has  presented  analytical  ideas  and  calculations  so  far 
accomplished.  It  has  also  presented  preliminary  ideas  for  analog  pressure 
distribution  calculation  and  display.  It  is  felt  that  both  analogs  and 
mathenmtical  analysis  have  much  to  offer  in  gaining  an  understanding  of  various 
features  of  the  near  sound  field  of  transducers  and  arrays.  The  analog  approach 
should  allow  a  rapid  semi-quantitative  assessment  of  the  value  of  various 
design  ideas,  with  no  limitations  on  transducer  shapes  and  distribution. 

The  analytic  approach  will  continue  to  give  the  necessary  deep  insight  into 
the  fine  details  of  array  geoMtries  that  are  sufficiently  simple  for  their 
representation  in  convenient  mathesmtical  expressions.  The  developsient  of 
analytic  approximation  techniques  should  be  continued,  the  objective  being 
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useful  engineering  formulae.  Another  important  phase  of  future  study  should 
be  a  study  of  the  applicability  of  x-ray  crystallographic  mathematics  for 
sound  array  problems.  This  should  teach  us  tricks  and  limitations. 
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1.  Analyli 

In  the  fourth  section  of  Pert  I  we  were  concerned  with  “.he 
evaluation  of 

(1)  I  -  4  f  ,  ,  y - 

where  previously  (“fc.  -  was  defined  as 

M 

(?)  S,(4'X^- '‘Vj  <5;/,r^)) 

Substitution  of  (2)  In  (l)  gives 
•• 

(3)1  -  :  j f  dwJw. 


Let  %  s  itxr  yfJ^nr^TTT^ 

aw  »■♦■>»•■  4.v»- 


xix  *  (»Tr)*’T>  Jv 


JiiL 

Om)‘V 


X*-- («»)•■( I)*-) 

»^-UTr)^CA^4.u-) 

CaT)*- 


■ji  g.  -V^ 

4*r 


(M 


T  »  Ai/iry  ^  - 

•  1  /cr  -i—  r«i.\».ircr-=“ 


ar  Uir)»-/Ri::;;  't 


aw- 
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Reference:  GrSbner  Wolfgang  und  Hofreiter  Nikolaus, 

Zwelter  Tell,  Bestlmmte  Integrale",  Sprlnger-Verlag, 

•e 


{ _ hh-i* _ * 

A 


"Integraltafel- 
1950,  p.  129,  No.  77c. 

>  o . 


We  have 


Let 


then  ^  ^ 


t  -  -1. 


A 


^  s  I  vrfV>ev>  *  • 
$=1^  wl>ew 


Then 


A  V  Sis  at 

Substituting  this  now  Into  Equation  (U) 

«• 

I  •  I  \  5l»i  (V«|)  ^  y  )  ‘ 

From  Oberhettlnger  p.  16^  ^  1^ 

-It  A  VSu 


0 

* 


This  agrees  with  the  results  that  were  obtained  by  merely  treating 
as  the  Dirac  delta  function  between  the  limits  of  and  •» . 


II.  Sumnary  of  Formulas 

\  rr  /■  p  \ 

1.  <t'»(5A;Ox[K(V>)'  •“'■r 

2.  •*«<*» 

-•e 
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3.  »■  “  1  * 

nTs*-**- 


-airi 


wkcrc 


4. 


5. 


6. 


7. 


8. 


9. 

10. 


•o 

(a*' >,♦*"->»;  43w(k.|^)w*’ (x'yjJnJ^. 

«»•> 

wUc 

S3CX-X')  «  Sin  (**^)  5;n  (»'^) 

(For  a  circular  piston) 


where 


4)  -  ..’.L, 

e 

S^cCair  )*•(>%>**’) 

s  =  Uir)’-CA%>^Vii) 

a 

V,  a  WCX,^)  wIlCM  O*’ 
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11.  (For  a  circular  piaton) 

i)  «  «.W.  \  —  3;  (r^)  T  (ra)  Jy* 

where  '*4=  wCx,^)  wVicm  ** 

/*  *  Vr‘.  -fc* 

V  « 

12.  (For  a  rectangular  piaton) 

C.sC»itAS)  •  H*lLVia*--CAiii)*- JciA 

_ A  w  A 


111.  Crapha  and  Preaaure  Contours 

1.  Figures  7*19  ahow  v(v)  calculated  for  the  thirteen  field  points 
for  the  sample  nlne-plston  array  used  In  the  study  of  the  accuracy  of  the 

generator.  The  horlsontal  scale  Is  measured  In  units  of  a,  the 
radius  of  the  pistons,  measured  with  respect  to  the  particular  field  point 
as  origin. 

2.  Figure  20  shows 

with  the  Parke  and  Moran  simulation  of  analog  computation.  Ihe  horlsontal 
scale  Is  measured  In  units  of  field  point  separation  CL  . 


comparing  the  Sherman  and  Kass  calculation 

fT  ^  ^ 
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